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Abstract. By using a time slicing procedure, we represent the solution operator of a second-order parabolic 
pseudodiflferential equation on R" as an infinite product of zero-order pseudodifferential operators. A similar 
representation formula is proven for parabolic differential equations on a compact Rieinannian manifold. Each 
operator in the multi-product is given by a simple explicit Ansatz. The proof is based on an effective use of the 
Weyl calculus and the Fefferman-Phong inequality. 
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1 . Introduction and notation 

We begin with recalling standard notation for the calculus of pseudodifferential operators (i/'DOs). 
Throughout the article, we shall most often use spaces of global symbols; a function a € "^^"(R" x R'^) 
is in 5"'(R" X R'') if for all multi-indices a, /3 there exists Ca/s > such that 

(1.1) |5'J5^fl(x,f)|<C„;;(^r-^l, .^eR", ^eR", {^} :^ (I + 

We write S'" = ^"'(R" x R"). i/'DOs of order m, in Weyl quantization, are formally given by (see [Hor79] 
or[H6r85, Chapter 18.5]) 

Op"(a) u{x) = a''ix, D,,)u(x) = (2ny" JJ e'<"'""'^>a((x + y)/2, ^) u{y) dy d^, u e ^'(R"). 

We denote by 1'"'(R"), or simply by the space of such i/'DOs of order m. 

We consider a second-order tfrDO defined by the Weyl quantization of q{x,^). Assuming uniform ellip- 
ticity and positivity for q{x, ^), we study the following parabolic Cauchy problem 

(1.2) d,u + q-\t,x,D,)u = 0, < f < r, 

(1.3) u |,=o = Mo, 

for Mo in L^(R") or in some Sobolev space. The solution operator of this Cauchy problem is denoted by 
U{t' , t), < t < f < T . Here, we are interested in providing a representation of t/(f', f) in the form of a 
multi-product of i/'DOs. 

Such a representation is motivated by the results of the second author in the case of hyperbohc equations 
[Le 06, Le 07]. If the symbol q is only a function of ^, the solution of (1.2)-(1.3) is simply given by means 
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of Fourier transformations as 

Following [Le 06], we then hope to have a good approximation of u{t, x), for small t, in the case where q 
depends also on both t and x: 

u{t, X) - (In)-" jj e'<.v-v,f>^-r,(0,(.+>)/2,f) ^ ^.v^^(^^ 

where p(,"f)(x,^) :- £-('"-'')?('■', r/)^ < f' < r" < T, which is in 5'". The infinitesimal approximation we 
introduce is thus of pseudodifi'erential nature (as opposed to Fourier integral operators in the hyperbolic 
case [Le 06]). 

With such an infinitesimal operator, by iterations, we are then led to introducing the following multi- 
product of i/rDOs to approximate the solution operator U{t', t) of the Cauchy problem (1.2)-(1.3): 

Pla){x,D,) ifO<f<f<", 
1 

Pl,4x, Z),)f[/7-,,. ,,,_„/x, D,) if f« < t < f(^-i). 

where ^ = f*", . . . , f<^') is a subdivision of [0, T] with = f*"' < < ■ ■ • < = T. It should be 
noted that in general pj'^',, ^,^{x, D^), t' < t", does not have semi-group properties. 

In [Le 06], for the hyperbolic case, the standard quantization is used for the equation and for the ap- 
proximation Ansatz. However, in the present parabolic case this approach fails (see Remark 2.7 below). 
Instead, the choice of Weyl quantization yields convergence results of the Ansatz "W^^,, comparable to those 
in [Le 06, Le 07]. The convergence of 'Wvp / to the solution operator U{t, 0) is shown in operator norm with 
an estimate of the convergence rate depending on the (Holder) regularity of q{t, x, ^) w.r.t. the evolution 
parameter t. See Theorem 3.8 in Section 3 below for a precise statement. 

Such a convergence result thus yields a representation of the solution operator of the Cauchy problem 
(1.2)-(1.3) by an infinite multi -product of i/rDOs. The result relies (i) on the proof of the stability of 
the multi-product "Wpj as N - |''^| grows to oo (Proposition 3.1) and (ii) on a consistency estimate that 
measures the infinitesimal error made by replacing U{t', f) by pj'^', ^^{x, D^) (Proposition 3.6). The stability in 
fact follows from a sharp Sobolev-norm estimate for p'^ ^^{x, D^) (see Theorem 2.2): for s e R, there exists 
C > such that 

(1.4) iip;;',,,)(x, DM(H-,.Hn < 1 + c(t' - 1). 

The Fefferman-Phong inequality plays an important role here. 

The representation of the solution operator by multi-products of i/^DO follows from the exact conver- 
gence of the Ansatz "Ws-^j in some operator norm. We emphasize that the convergence we obtain is not 
up to a regularizing operator A further interesting aspect of this result is that each constituting operator 
of the multi-product is given explicitly. With such a product representation, we have in mind the devel- 
opment of numerical schemes for practical applications. Once the problem is discretized in space, the use 
of fast Fourier transformations (FFT) can yield numerical methods with low computational complexity, 
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with possibly microlocal approximations of the symbols in question as is sometimes done in the case of 
hyperbolic equations (see for instance [dHLWOO, LdHOla, LdHOlb, LdH03]). We also anticipate that our 
representation procedure can be used in theoretical purposes. 

As described above, the first part of this article is devoted to the parabolic Cauchy problem on R" and 
to the study of the properties of the approximation Ansatz "W^^j. In the second part, we shall consider 
a parabolic problem on a compact Riemannian manifold without boundaries. In this case, the operator 
q'^ix, Djc) is only considered of differential type, in particular for its full symbol to be known exactly. In 
each local chart we can define an infinitesimal approximation of the solution operator as is done in R" and 
we combine these local i/rDOs together with the help of a partition of unity. This yields a counterpart of 
/jj^, j^{x,Dx) for the manifold case (see Section 4.1) denoted by P{t'j). In fact, the sharp estimate (1.4) still 
holds in L? {s - 0) for this infinitesimal approximation (Theorem 4.4). The proof of a consistency estimate 
(Proposition 4.6) requires the analysis of the effect of changes of variables for Weyl symbols of the form of 
-t)q^ -pjjg choice we have made for the definition of P(t',t) is invariant through such changes of variables 
up to a first-order precision w.r.t. the small parameter h - t' - 1, which is compatible with the kind of results 
we are aiming at. With stability and consistency at hand, the convergence result then follows as in the case 
ofR". 

In the manifold case, the constituting i/rDOs of the multi-product are given explicitly in each local chart. 
We observe moreover that the computation of the action of these local operators can be essentialy performed 
as in the case of R", which is appealing for practical implementations. 

Another approach to representation of the solution operator U{t', t) can be found in the work of C. Iwasaki 
(see [Tsu74, Iwa77, Iwa84]). Her work encompasses the case of degenerate parabolic operators, utilizes 
multi-product of i/^DOs and analyses the symbol of the resulting operator, using the work of Kumano-go 
[Kg81]. However, the symbol of the solution operator U{t',t) is finally obtained by solving a Volterra 
equation. Such integral equations also appear in related works on the solution operator of paraboUc equa- 
tions (see e.g. [Gre71, ST84]). The alternative method we present here will be more suitable for applica- 
tions because of the explicit aspect of the representations. The step of the integral equation in the above 
works makes the representation formula less explicit. However, the reader will note that the technique we 
use in our approach here do not apply to the case of degenerate parabolic equations like those treated in 
[Tsu74, Iwa77, Iwa84]. The question of the extension of the convergence and representation results we 
present here to the case of degenerate parabolic equations appears to us an interesting question. 

Let us further recall some standard notions. We denote by cr(., .) the symplectic 2-form on the vector 
space r*(R"): 

(1.5) cr((x,a(y,r]))^{^,y)-{r],x}, 
and we denote by {/, g] the Poisson bracket of two functions, i.e. 

n 
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We shall use the notation #" to denote the composition of symbols in Weyl quantization,i.e., a" (x, Dx) o 
b"'(x, D,) = {a #"bY{x, D,). The following result is classical. 

Proposition 1.1. Let a e S"\ b e S'"'. Then a #"b e 5'"+'"' and 
(1.6) 

k 



(a r'b){x, O^Y- {^(riiDx, ), (D,„ D,))) a(x, ^)biy, rj) I 
1 ^ 



® 
where ^, f , r, ^) = 2«t - ^, z) - - ^, f )). 

The result of Proposition 1.1 is to be understood in the sense of oscillatory integrals (see e.g. [H6r90, 
Chapter 7.8], [AG91], [GS94] or [Kg81]). For the sake of concision we have introduced 



for« > 3, n € N. 

® 



For the exposition to be self contained, we prove Proposition 1 . 1 in Appendix A. 

We sometimes use the notion of multiple symbols. A function a(x,^,y, rj) e C°°{W^ x R''' x R''- x R''-) 
is in 5"''"''(R'" x R'" x R''^ x R''^), if for all multi-indices q;i,/3i, a2,y82, there exists C^',^^ > such that 

x e R«', y e W\ ^ e R'", e R''= (see for instance [KgSl, Chapter 2]). 

For s e R. We set E^''> := {D_,y = Op((^)'), which reaHzes an isometry from H''{R") onto H'-'iR") for 
any r e R. We denote by (., .) and ||.|| the inner product and the norm of L^(R"), respectively and || . \\fj' for 
the norm on H%R"), s eR. For two Hilbert spaces K and L, we use || . \\(k,l) to denote the norm in X.{K, L), 
the set of bounded operators from K into L. 

Our basic strategy is to obtain a bound for i/rDOs involving a small parameter h > 0. In the following, 
we say that an inequality holds uniformly in h if it is the case when h varies in [0, /Jmax] for some /i^ax > 0. 
In the sequel, C will denote a generic constant independent of h, whose value may change from line to line. 
The semi-norms 

(1.8) Pa^a)-- sup {^}-'"^^)d1cfMx,0\ 

(A-,f)eR"xR'' 

endow a Frechet space structure to S'"(R" x R''). In the case of a symbol a/, that depends on the parameter 
h we shall say that a;, is in 5™^ uniformly in h if for all a, (3 the semi-norm Pajsiai,) is uniformly bounded 
in h. Similarly, we shall say that an operator A is in ^P"" uniformly in /; if its (Weyl) symbol is itself in 5'"' 
uniformly in h. 

The outline of the article is as follows. Sections 2 and 3 are devoted to the multi-product representation of 
solutions on R" . In Section 2 we prove the sharp Sobolev norm estimate ( 1 .4), which leads in Section 3 to the 
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Stability of the multi -product representation. We then prove convergence of the multi-product representation 
in Section 3. Some of the results of these two sections make use of composition-like formulae, whose proofs 
are provided in Appendix A. In Section 4, we address the multi-product representation of solutions of a 
second-order differential parabohc problem on a compact Riemannian manifold. As in the previous sections 
we prove stability (in the case) through a sharp operator norm estimate and we prove convergence of 
the multi-product representation. The convergence proof requires an analysis of the effect of a change of 
variables on symbols of the form e from one local chart to another, which we present in Appendix B. 

2. A SHARP //' BOUND 

We first make precise the assumption on the symbol q{x, ^) mentioned in the introduction. 

Assumption 2.1. The symbol q is of the form q — q2 + q\, where qj e S^, j — 1,2, q2(x,^) is real-valued 
and for some C > we have 

q2{x,^)>C\^f, x€R", ^ eli", \^\ sufficiently large. 

Consequently, for some C > 0, we have 
(2.1) q2ix, ^) + R&qy(x, ^) > C\^\-, x E R", ^ e R", |^| sufficiently large, say |^| > i? > 0. 
As is stated in the introduction, our main aim is to deal with the operator p''i^{x, D^) where 

p„(x,^)^e->'''^^-^. 

It is well-known that the i^DO p]l(x, D^) is uniformly i/* -bounded in h, i e R. Actually we have the 
following sharper estimate. 

Theorem 2.2. Let i e R. There exists a constant C > such that 

\\pl{x,DMiH',Hn<^+Ch, 

holds for all h > 0. 

To prove Theorem 2. 1 we shall need some preliminary results. 
Lemma 2.3. (i) Let I > and r e S'. Then h'^^r pi, is in uniformly in h. 

(ii) Let a and P be multi-indices such that \a + j3\ > 1. Then, for any < m < 1, we have d"cf^ph - h'"p"i'"^, 
where p"^"^ is in S uniformly in h. 

Proof. We have 

(/i<^2))./g-/>Re9(x,f) < f.,^ xeR", ^eR", h>0, 
for all y e N by (2.1), hence 

h"\ix, OPhix, ^\<C, xeR", ^ e R", h> 0. 
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For multi-indices a and /3 we observe that h'l^d"dP^{r{x,^)ph{x,^)) is a linear combination of terms of the 
form 

h'^^>l\d'',°S^'r){x,^){d",^dP^q){x,^) ■ ■■{d",'dP^q){x,^)pH{x,^), 

for A: > 0, ffo + H ha^ - a, Pq + [5\ -\ v Pk - [i and the absolute value of this term can be estimated 

by 

//■+'/2(^>'+2*-^I|/?a(x,^)| <C(^r^l, jceR", ^eR", h>Q, 

by (2.1), which concludes the proof of (i). For \a + 13\ > 1, d"d^ph(x,^) is a linear combination of terms of 
the form 

for ^ > 1, ffi + ■ • ■ + fft - a, Pi + ■ ■ ■ + Pk - P, which can be rewritten as li"'Ai,{x, where 

i„(x,^) = (d:'d^^'q)(x,o ■ ■■(d:'d^^'q){x,mf""^'{m^f Phix,o- 

Since (5"' 5^'^) ■ ■ ■ (d"/d'^'q) e S^'^-^K we see that A,, e S'"'-^^ uniformly in h by using (i). ■ 
From Weyl Calculus and the previous lemma we have the following composition results for the symbol 

Ph- 

Proposition 2.4. Let r/, be bounded in S', I uniformly in h. We then have 

(2.2) r;, #'"ph = rhPh + h^-/°'' = r/,/?;, + /z/lj^'^ = r/,/?;, + ^{r/,,/?/,} + hAf\ 

(2.3) Ph #'V/, = rhpi, + hipf^ = r/,/9/, + /i^u}/' = ri,ph + Yi^Ph, n,] + hjjif\ 

where A^i^\ ~A^^\ and fif^ are in S' uniformly in h and /ij^'' and /ij^'' are in 5'^' uniformly in h. 

To ease the reading of the article, the proof of Proposition 2.4 has been placed in Appendix A. We apply 
the result of Proposition 2.4 to prove the following lemma. 

Lemma 2.5. We have <^)^' #">/, - <^)^' \ph\^ - hkk with k/, in S uniformly in h. 

Proof. By Proposition 2.4 we have 

1 1 " 

i^f #''Ph = i^f Ph + [{^f\Ph] + hAi,H = i^f' Ph + Y- ^-^.-P" + ''-^1.'" 

with Ai j, in S^" uniformly in h. We then obtain 

pT, #"' (f)'-' #>„ = pT, r {{^f' PH + ^ J(%(^)'0 5,^;./,) + /z^2,A, 
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with A2,h in uniformly in h. By Proposition 2.4 we have 

W^ (<^)'' Ph) = (O^'lPhf + Yi <^>''/'"} + ''^ 

1 " 

with /l3_/, in 5^* uniformly in /j. We also have 

with A41, in 5^', uniformly in h, by Proposition 2.4. We have thus obtained 
(2.4) #"■ <^)-' = {^-'IPhf + ^4" + 5J Z + 

with /I5 J, in 5 "^^ uniformly in /j, and with 

n 



and 



We introduce a := ^2 + Re <?i and yS := Im^i. We have 
(2.5) 4" - 2//!V/.l'(^>'' 

and 

if] = -2//z|p„|2(5,/) %<^)2'. 

Since a e and /? € 5' we then have {a,/3} e S^. From Lemma 2.3, we thus obtain /J^'' = with 
in 5^' uniformly in h. We also have that \ - hkf\ with kf \ in 5^' uniformly in h, which from (2.4) 
concludes the proof. ■ 

We shall also need the following lemma. 

Lemma 2.6. We have {^Y#'' \ph?#'' {^Y - <^)->aP = hkh with kh inS^' uniformly in h. 

Proof. We set p/, = {^Y Ip/.p #»' {^Y- From Weyl calculus we have 

Ph{x, = 7r-2"J e'2fe.'.f'^'f><^)'<T)^ + z + f, f) fifz dt dT 

® 
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where t, t, ^) - 2({z, t - ^) - {t,( - ^)). Arguing as in the proof of Proposition 1 . 1 in Appendix A we 
write 



Ph{x, ^)-{4f'\pH\\x, a = 7r-2«J e'^^-''^'''^\0\TY {\ph\\x + z + t,^)- \ph\\x, 0) dz d( dt dr 

® 

(2.6) = '-n-^" 2 // '("5-. - di)({0\ry) + r{z + t), ^) dr dz d( dt dr, 

■'■=' ® 

by a first-order Taylor formula and integrations by parts w.rt. ^ and r. Observing that we have 
° = ? Z^^^' " '5^,)«^)'<T>0 {d,j\ph\^){x,0 |^_^_^ 

i 

we can proceed as in the proof of Proposition 1 . 1 (integration by parts w.rt. r in (2.6) and further integrations 
by parts w.rt. C, and r) and conclude after noting that \ph^ satisfies the properties listed in Lemma 2.3 like 
Ph- ■ 



Remark 2.7. Note that the use of the Weyl quantization is crucial in the proofs of Lemmata 2.5 and 2.6. The 
use of the standard (left) quantization would only yield a result of the formp^#(^)^* #/?/, = h^-kh 

with kh in S uniformly in Such a result would yield a h'- term in the statement of Theorem 2.2 and the 
subsequent analysis would not carry through. 



We now define the symbol v/,(x, ^) = idMiifl^ for h > 0, and prove the following lemma. 
Lemma 2.8. The symbol v/, is in S ^ uniformly in h. 

Proof. We write v/,(x,^) = 2 Req{x,^) g-2'/iRe?(.i-,f) -pj^g integrand is in 5'^ uniformly in r e [0, 1] and 
h by Lemma 2.3 and Req e S^. ■ 

Lemma 2.9. The symbol ph is such that 

{{\ph\^nx, D,)u, m) < (1 + Ch)\\uf, u e L\R"), 
for C >0, uniformly in h > 0. 

Proof Let x e %?°{R"), < ;r < 1- be such that x(0 = 1 if 1^1 ^ ^- Then we write 

\ph\Hx,0 = e-2/,Re<,(x,f) ^ ^-2Mf)Re<?(..f)^-2Ml-^(f))Re<,(x,f) ^ (1 + /,^(_,, ^))e-2''(>-^«»R«?(-.?), 

where yL/(x,^) = -2;t(^) Re ^(jc, ^) j^' g-2rhx(m<=<iU4)j^^ Prom [Hor85, 18.L10],the symbol is 
in 5 ° uniformly in r and h; hence the symbol yu(x, ^) is in 5 uniformly in h. From [Hor85, Theorem 18.6.3] 
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and Lemma 2.3, it is thus sufficient to prove the result with ph(_x,^) replaced by pi,{x,^) = e 
We set 

for h > 0. From (the proof of) Lemma 2.8 we find that v/,(x, ^) is in 5^ uniformly in h. Since 1 - xiO - 
if 1^1 < 1? and Ktq{x,^) > if |^| > we observe that Vh{x,^) > 0. Then the Fefferman-Phong inequahty 
reads ([FP78], [H6r85, Corollary 18.6.11]) 

[vlix, D,)u, u) > -C\\u\\l„ u 6 l\R"), 

for some non-negative constant C that can be chosen uniformly in h. This yields 

M\l2 - {{\phfnx,D,)u, u) > -Ch\\u\\i„ u 6 L-(R"), 

which concludes the proof. ■ 

We are now ready to prove Theorem 2.2. 
Proof of Theorem 2.2. We use the following commutative diagram, 

W > W 

l2 !i , l2 

and prove that the operator T/, satisfies ||r/,||(£2^2) < 1 + Ch. 
The Weyl symbol of T^* o T/, is given by 

(Th = <^)-' r pj; r {^f' r'p„ r <^>-. 

By Lemmata 2.5 and 2.6, we have cr/, = \pi^'^ + hk/,, with k/, in 5" uniformly in h. We note that kh{x,^) is 
real valued. To estimate the operator norm of T/, we write 

\\Thu\\- - (jl o TijU, m) = {i\ph\~T'ix, Dx)u, ii) + h {klix, Dj,)u, uj 
<{(\pi,\^r'(x,D,)u,u) + Ch\\uf, 

for C > [Hor85, Theorem 18.6.3]. The result of Theorem 2.2 thus follows from Lemma 2.9. ■ 
Let m{x) be a smooth function that satisfies 

(2.7) < OTtnin < m(x) < mmax < <», 

along with all its derivatives. With such a function m, we define the following norm on L^(R") 

II/IIl^(R..,„,..) = / f(x)'n{x)dx, 

R" 

which is equivalent to the classical L? norm. We shall need the following result in Section 4. 
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Proposition 2.10. There exists a constant C > such that 

\\pI{x, £>.v)m|Il2(R",,„^/.i-) < (1 + Ch) ||M||^2(R„_,„rf^), u € L^(R"), 

holds for all h > 0. 

Proof. We follow the proof of Lemma 2.9 and use ph{x, ^) in place of ph{x, We then set 

~y,(x,^):^'!^{l-\p,f(x,0). 
Then Vh{x, ^) > is in 5 ^ uniformly in h. The Fefferman-Phong inequality yields 
{y';;(x, D,)u, u) > -C\\u\^, > -C'||m||2,(^„ u e L^CR"). 

This yields 

II«IIl^(R",„,..) - [iMPh^nx, D^)u, u) > -C'/,||«|l^,(^„ u e 

By Lemma 2.11 just below, we have 

[{m\ph?-T\x, D,)u, u) = iji^\x, D,) omo {pl{x, D,)u, u) + hiX^ix, D,)u, u), 

with A{x,^) in 5" uniformly in h and where ni stands for the associated multiplication operator here. We 
conclude since pi^{x,Dx) - (p''l{_x,D_^))* . ■ 

Lemma 2.11. Let f e "^""(R") be bounded along with all its derivatives. We have 

'p'h#''f#"-'ph - f\ph\' = h'^h, 

with Ah in uniformly in h. 

Proof. From Proposition 2.4 we have 

1 " 

/#>/, = fph - Yi ^^^-'if^ + '^'^1.'" 

with Aij, in 5" uniformly in h. By Proposition 2.4 we also have 

Ph#''(fph) = flphf + -[phjph] + hAlJr 
2i 

= f\Ph?- + 5^ X {(^?jPh)(.d-^jf)Ph + id^jPh)fidjcjPh) - id^,:ph)fidi;jph)) + hA2,h, 
M 

with /I2,;, in 5° uniformly in h, and 

'Ph#" {(dxjf)d£jPh) = 'Ph(9.xjf)d^jPh +hfij^i,, j ^\,...,n, 

with fij i, in 5° uniformly in h. It follows that 

1 " 1 
Ph#"f#''ph ^ flphf - ^ -phd^^'pD + — f {Ph, Ph] + hA4,h, 

with /I4,;, in 5" uniformly in h. With the notation of the proof of Lemma 2.5 (see expression (2.5)) we have 

{pruPh]^2ih\i,f {a,p]^hk'-^\ 
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with A:^'^ in 5 " uniformly in h by Lemma 2.3, since yS is in 5 ' and a e 5 ^. We also have 

n 

with fcjj^' in 5" uniformly in /i by Lemma 2.3. ■ 

3. Multi-product representation: stability and convergence 
We are interested in a representation of the solution operator for the following parabolic Cauchy problem 

(3.1) d,u + q"\t, X, D,)u = 0, < f < T, 

(3.2) M |,=o = «o e H\R"). 

Here the symbol q{t,x,^) is assumed to satisfy Assumption 2.1 uniformly w.rt. the evolution parameter t 
and to remain in a bounded domain in 5 ^ as f varies. We then note that the result of the previous section 
remains valid in this case, i.e., the constant C obtained in Theorem 2.2 is uniform w.r.t. t. We denote by 
U{t', t) the solution operator to the evolution problem (3.1). 

Following [Le 06], we introduce the following approximation of U{t, 0). With ^ - {f^°\ f^'\ . . . , f'^^), a 
subdivision of [0, J] with = f'°* < f^'' < ■ ■ • < f^'^' - T, we define the following multi-product 

^(,,0) if < f < f">, 

1 

P{,,t"")Y]P{tV\,v-n) if f<^'' < t < f<*+». 



(3.3) nV.^,, :-- 



where P(t\i) is the ifrDO with Weyl symbol p(,',t) given by := e"'' -t)q(t,x,£) j^^j. ^> ^ 
Pi,', Ax) = ,/x,D,)y(x) = {In)-" jj e'<-v/>g-(''-0?(a.v-+y)/2,f,^(^) 
We shall prove the convergence of to U (f , 0) in some operator norms as well as its strong convergence. 
3.1. Stability. As a consequence of the estimate proven in Theorem 2.2 we have the following proposition. 
Proposition 3.1. Let s e R. There exists K >0 such that for every subdivision of[0, T], we have 

Vf€[o,r], |inv-i,,iU-.ff-) ^ e'^^- 

Proof. By Theorem 2.2, there exists C > such that we have \\P(,',,)\\^fj, < 1 +C(t'-t) for all f, t e [0, T], 
t' > t; we then obtain 

N-l 

\\'W%\r\\^H',H')^Y]i^+C{t^'^'^ -t'-'^)). 

i=0 

Setting = ln(nZo'(l + C(f('+" - f^''))), we then have t/.;, < Xf=o C(f^''+'' - f^'*) = CT. We thus obtain 
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3.2. Convergence. To obtain a convergence result we shall need the following assumption on the regularity 
of the symbol q{t, x, ^) w.r.t. the evolution parameter t. 

Assumption 3.2. The symbol q(t, x, ^) is in '^°'«([0,r],5-(R" X R«)), i.e.. Holder continuous w.r.t. t with 
values in S^, in the sense that, for some < ff < 1, 

q{t', X, ^) - q{t, X, ^) = (f' - 0" §(f', t, x,^), Q<t<t' < T, 

with q(t' , t, X, ^) in S ^ uniformly in t' and t. 

We now give some regularity properties for the approximation Ansatz "Ws-^j we have introduced. 

Lemma 3.3. Let i e R and t" , t e [0, T], with t < t". The map t' i-» P(t'fy for t' e [t, t"], is Lipschitz 
continuous with values in J2{H^{R"), H^^^{R")). More precisely there exists C > such that, for all v e 
H\R") andt^^\t^^'> e [t,t"]. 

Proof. We simply write 

,(2) 

- P(,(„,o)(v)W = -(In)-" j j'J e'<-v.f>e-(''-'MU.+j)/2,f) q^f^ + yy2, v{y) dy d^ dt' . 

/(') 

We thus obtain a if/DO whose Weyl symbol is in 5^ uniformly w.rt. f*-* and f*'' and conclude with Theorem 
18.6.3 in [Hor85]. ■ 

Lemma 3.4. Let s e R, t" , t e [0, T], with t < t", and let v e H\R"). Then the map t' i-> f(,',,)(v) is in 
'^^\[t, t"],H'(R")) n '^'([f, t"],H'-HR")). 

Proof Let f*'> e [f, t"] and let e > 0. Choose vi 6 H'+^R") such that ||v - viH^, < s. Then for f*^* e [f, t"] 

||P(,a,.,)(y) - P(,m,,p)\\^, < ||^'(,(2),,)(v - vi)||^,, + ||f'(,(i),,)(v - vi)||^, + ||^'(,a,.,)(vi) - P(,<M,oivi)\\^, 
(3.4) < 2(1 + C(t" - t))e + C |r(2) _ ||v,||h-2. 

The continuity of the map follows. Differentiating P{i\i){v) w.r.t. f , we can prove that the resulting map t' i-» 
di'P(i'j)iv) is Lipschitz continuous with values in X.{H''^^(R"), H^'^(R")) following the proof of Lemma 3.3: 
there exists C > such that for all w e H'^^(R") 

||(5,P(,(2,,,) - d,,P(,o>.,))M\\^,-2 < Cy^^ - ff'^l ||w||h-2. 

Here dt'P{tV)j) means We also see that the map v i-> dt'P(t',t)(v) is continuous from H\R") 

into //^"^(R") with bounded continuity module: with v € H\R"), we make a similar choice as above for 
vi e H^^^(R") and obtain an estimate for 

||5,'^(,(2.,o(v)-5,P(,u),,)(v)||^,_, 

of the same form as in (3.4). ■ 
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Gathering the results of the previous lemmata we obtain the following regularity result for the Ansatz 

Proposition 3.5. Let 5 e R, let mq e i/^(R")- Then the map "W^^jiuo) is in "^•^([O, T], H'iW)) andpiecewise 

The following energy estimate holds for a function /(?) that is in ^°([0, T], //'(R")) and piecewise 
'g''([0, T],H''\W)) (by adapting the proof the energy estimate in Section 6.5 in [CP82]): 

T T 

(3.5) 11/(011^.-,+ J Wmwl, dt' < C ||/(0)||^,_, + J Udr+q''(t',x,D,))f(t')\\l,_,dt' , 


for all t 6 [0, T]. Once apphed to (t/(f, 0) - nV>p,,)(Mo) with uo £ //'(R") we obtain 

T T 

\\(U(t, 0) - nV.iv)(Mo)llH.,-. + J ll(f^(f> 0) - nV.iv)(Mo)llH. dt'<cj \\ {dr + q'\t',x, D,)) nV.iV'("o)llH,,-: 



= C 2 j II (5,, + q'\t', X, DJ) P(r,,o))nV.^},,o,("o)llH,,-2 fl^f' 

(3.6) < C 2 j II (a, + q'^if, X, £),)) P(r,,<i))ll(H., H.-) "'f' ^''^INolli,, 

where we have used the stability result of Proposition 3.1. It remains to estimate the Sobolev operator 
norm of (5,' + q'''{t' , x, Dx)) P(t',i), for t' > t, which can be understood as estimating the consistency of the 
proposed approximation Ansatz. This is the object of the following proposition. 

Proposition 3.6. Let s eR. There exists C > such that \\{d,' + ^"'(f', x, D j.))P(,/_,)||(^,^,,-2) < C(f' - f)" for 
0<t<t'<T. 

Proof. We have 

(3.7) ^,P^,,,^uix) = -iln)-" jj e'<^-^''f>^(f, {x + y)ll, ^)e-(''-')?(',(x+y)/2,f) 
and thus the operator {df + q^^\t' , x, D J)P(,' ,) admits 

o-(t',t) = q{t', ., .) #>(,',o - q{t, ., .)p(t',t) 
for its Weyl symbol. Since by Assumption 3.2 we have 

q{t', X, ^) - q{t, X, ^) = it' - tf q{t', t, x,^), < t < t' < T, 
with q(t' , t, X, ^) in 5 ^ uniformly in t' and t. We can thus conclude with the following lemma. ■ 

Lemma 3.7. We have q #"'ph — qpi, — hA/, with Ah in S'^ uniformly in h. 
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Proof. As in Section 2, we ignore the evolution parameter / in the notation. The result is however uniform 
w.r.t. t. By Proposition 2.4 we have 

q #"'ph = qph + ^.lq,Ph] + h^^h, 
2i 

with Ai, in 5^ uniformly in h. We note however that {q, pi,} = 0. ■ 
The result of Proposition 3.6 and estimate (3.6) yield 

(3.8) Wmt, 0) - nV-i,,,)(Mo)llH..-. + (J \mt, 0) - ^5i,r)(Mo)|li> t/fj' < CTe" \\uo\\h', 



where As^* = maxo<y<Ar-i(f''^' - f'-'')- This error estimate implies the following convergence results which 
provides a representation of U{t, 0) by an infinite multi-product of t//DOs: U{t, 0) = hm^.^^o ^%!,f- We now 
state our main theorem. 

Theorem 3.8. Assume that q{t,x,^) satisfies Assumptions 2.1 and 3.2. Then the approximation Ansatz 
•Wi^i^f converges to the solution operator U(t, 0) of the Cauchy problem (3. l)-(3.2) in JUiH'iW), //'"'^'^(R")) 
uniformly w.r.t. t as A-|i = maxo<)<A'-i(f''^' - f^-'') goes to with a convergence rate of order a(\ — r): 

\YW.^,t - U(t, 0)||(„.,H.,-..o < CA;*'-'-\ te[0,T], < r < 1 . 

The operator "W^^j also converges to U(t, 0) in L^{Q, T, X.(H\R"), H\R"))) with a convergence rate of 
order a: 

Furthermore 'W<^^i strongly converges to Uit, 0) in ^.{H^R"), //^(R")) uniformly w.r.t. t 6 [0, T]. 

Proof. The first two results are consequences of (3.8). The proof of the first result for r follows by 
interpolation between Sobolev Spaces [LM68]. 

Let Mo e H^{W) and let e > 0. For the strong convergence in H\R") we choose «i E //'^'(R") such that 
ll"o - miIIh' ^ s. We then write 

m^vAm) - t/(f,0)(M0)||H.. < ||W.?,,(mo - Ml)llff. + mvAui) - f/(f,0)(Ml)||H,> 

+ \\U{t,Q)(uo - Mi)||ff, < Ce + CA|||Mi||ff,.,, 

from the case r = of the first part of the theorem and from the stability of 'Wf^.j (Proposition 3.1). This 
last estimate is uniform w.rt. f E [0, T] and yields the result. ■ 

4. Multi-product representation on a compact manifold 

4. 1 . Notation and setting. We shall now consider the case of a parabolic equation on an «-dimensional 
compact "^^""-Riemannian manifold (At, g), where g is a smooth Riemannian metric. We let A be a second- 
order elhptic difl'erential operator on M whose principal part, A2, is given by the Laplace-Beltrami operator 
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on Al, which reads 

in local coordinates, where g = det(g,y). Other uniformly-elliptic operators can be considered by changing 
the metric. We choose here to focus on the differential case instead of the pseudodifferential case because 
the full symbol of the operator can then be completely defined on the manifold Ai. 

We allow the operator A to depend on an evolution parameter t. We shall thus assume that the metric is 
itself time-dependent, yet continuous w.rt. t, g = g(t, x), and satisfies 

(4.1) < c < g(f,jc) < C < oo, te[Q,T], x e M. 

For the L? norm on Ai, we shall use the metric g{0, x) as a reference metric. We set go(x) = g(0, x). We 

1/2 

then denote by dv the volume form which is given by dv — g^^ (x)dx in local coordinates. The L--inner 
product is then given by (u, w) = J^uwdv [Heb96]. 

Since we are going to consider an infinite product of i^DOs, a little attention should be paid to a finite 
atlas. We shall use an atlas = {Oi, ifrfjiei of Ai, \I\ < oo, with i/f; : 6i — > 6i, where 6i is a smooth bounded 
open subset of R". For / e J, we set 

Ji := € J; 0i n 0j + 0), J-f {I e Jf, j e J,}, 

which lists the neighboring charts and the "second"-neighboring charts for the chart (0,, (/f,). For technical 
reasons, we shall assume that there exists a coarser finite atlas = (0^, ^it)ke'K of Ai, *Pa: : 0^ — » 0a c R", 
such that for each chart (0,-, i/^,) 6 £/ there exists a chart (0a:(,), ^k(i)) £ such that 

i.e., 0jk,) contains all the "second"-neighbors of 0,. This is always possible by choosing the atlas suffi- 
ciently fine. We shall denote by fl,(f), ' e the Weyl symbol of A(t) in each local chart (fl,, if/j). 

We set {(fidiei as a family of real-valued functions defined on Ai such that the functions ((pf)ieT form 
a partition of unity subordinated to the open covering (6i)iei, i.e., 

supp((/j,) c di, < ifii < I, iel, and '^(pj - I. 

We denote 

Vi = (</'r')Vi = ° 

(21 

and similarly, for I ^ J^] , we shall set 

with ^k(i) as above, when there is no possible confusion on k{i). 
We set Q{t) as the elliptic operator on M defined through 

A{t) = ^V/oeCOo^Z- 
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The construction of Q can be done recursively: we write Q - Qi + Qi + Qo, with Qi a differential operator 
of order /, Z = 0, 1 , 2 and obtain 

Qi^A, Qi ^ -Yj[<Pi,Q2]o(pi, Qo = -Yj[(Pi,Qi]oipi. 

iel iel 

The recursion stops after two iterations since we consider differential operators here. 
In each local chart (0,-, xj/i), i e J, we denote by ^,(r, x, ^) the Weyl symbol of Q{t), i.e., 

Vm e %°°iei), Q{t)u = ri [q^XU X, D,){{ilf.'Yu)) , 

or equivalently 

V5 e qjit, X, D,,)u = {ilf.'r mWiU)) . 

The symbol ^,(f , x, ^) is uniquely defined since Q{t) is a differential operator. We also let quit, x, ^) be the 
Weyl symbol of Q{t) in the chart (0*, ^u), k e'K. From (4. 1) we then have 

Lemma 4.1. In each chart the symbol of Q(t) satisfies the properties of Assumption 2.1. 

We set 

Pui',f){x, ^) = iel, 0<t<t' <T, xedi, R". 

With these symbols in S ^\9j x R"), we define the following t/^DOs on M: 

(4.2) Pi,(t',t)U -.^ipiO if/* o p^,. ,^ix, Dt) o (i//V^y o ipi = (/r* o (ijpi o /?,■ (,, , Z).v) o ipi) o (l/r^')*, 

(4.3) := 2 p,. 

iel 

where (p, and ^, are understood here as multiplication operators. The operator P(t'.t) is the counterpart of 
the operator ^^{x,D^) introduced in Sections 2 and 3. We shall compose such operators in the form of 
a multi-product as is done in Section 3 to obtain a representation of the solution operator to the following 
well-posed parabolic Cauchy problem on M 

(4.4) d,u + Ait)u = 0, < f < r, 

(4.5) M |,=o = MO e H'iM). 

We denote by U{t', t) the solution operator of (4.4)-(4.5) and we define the multi-product operator ^>p,/ as 
in (3.3) for a subdivision ^^3 = f^'*, . . . , f^} of [0, T]: 

( Pm if < f < f('>, 

1 

P(rjm)Y]P(,iiK,ii-u-, if ff*^' < t < ff'^+l). 

We shall make the following regularity assumption on the operator A(f), which is equivalent to that made 
in Section 3 (Assumption 3.2). 



(4.6) ^9 
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Assumption 4.2. The symbol of A{t) is Holder continuous of order a, < a < 1, w.r.t. t with values in S ^: 
for each chart (0;, i/',) we have a,- 6 '^"•"'([O, T],S^(R" X R")), in the sense that, 

aiit', X, ^) - a,(r, X, ^) = (f' - tf ciiit', f, x,^), < t < t' < T, 

with ai{t' , t, X, ^) in uniformly in t' and t. Note that the same property then holds for the symbol ofA{t) 
in any chart. 

This property naturally translates to the symbols qi(t), i e I. 

Remark 4.3. The form we have chosen for the operator P(t',,) can be motivated at this point. First, a natural 
requirement is that = Id, which is achieved since Yjieifj - 1- Second, the consistency analysis of 
Proposition 3.6 gears towards having (d',P(fj) - A(f') o P(,>j))\^,_^ - 0, which is achieved here thanks to the 
form we have chosen for the differential operator Q{t). 

As in Section 3, we first need to address the stability of the multi-product. Here, we shall only consider 
the case. 

4.2. Stability. As in Section 2, we find a sharp estimate of the L^-norm of the operator over M. 
Theorem 4.4. There exists a constant C > such that 

\\P(i',i)\\(lHM).l^(M)) ^ 1 + C(t' - t), 

holds for all < t < t' < T. 

Therefore, as in Section 3, we obtain the following stability result for 'W-;!,,. 
Corollary 4.5. There exists K > such that for every subdivision ^ o/[0, T], we have 

W€[o,r], 

Proof of Theorem 4.4. We let u,w e L^{M). We have 

(P(,'-f)M- vv) = ^ I ^,M7i/'*(/?"(,,,)(jc,D,)((iA^')*((/;,M)))(iv= ^ I (ipiWi) p'-^,, ,^(x,D_,)(ipiUi) q\I^(_x) dx, 

where S,- (resp. vv,) is the pullback of u\b^ (resp. w\g.) by i/tt'. We now extend the symbol ^,(f, .) to R" x R" 
to obtain a symbol satisfying Assumption 2. 1 like its counterpart in Section 2. We still denote by ^,(f, .) this 
extended symbol. Then, by Proposition 2.10, for all / € J, there exists C, > such that 

where go is also extended from 0, to R", yet still preserving Property (2.7). With C - max,gj C, (recall that 
I is finite) we thus obtain 

I {P(f-,)u, w) I < (1 + c(f' - 0) Yj \\'f'>'^>\\L\0.,Q\r-d.^) ii'5^'"'iiL2w.gr -i^)- 
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Figure 1 . Change of variables bringing the analysis to the chart i&ka), ^k(i)) for the charts 
(ft, iffj) and {9j, ijjj), j e J'i, and their neighboring charts. 

A Cauchy-Schwarz inequality then yields 

I {P,r-,u, w) 1 < (1 + Cif - t)) ^i^^^_^.J (2 ^^^^^.^^J ■ 

iel iel 

Observing that 

iel " iel g_ iel ^ 

since 2,ej = 1, we find 

I {P(f-t)u. w) I < (1 + C(f' - 0) M\v-(M) M\n(M)^ 
which concludes the proof. ■ 

4.3. Consistency estimate. As in Section 3, Proposition 3.6, for the case of R", we shall now analyze 
the symbol of the operator {df + A(t'))P(,'j) and prove the following proposition that corresponds to a 
consistency estimate. 

Proposition 4.6. LetO <t <t' <T. We have 

(d, + A{t')) o = (f - f)"L(,,,), with L(,v) e "VHM), 
and for all 5 € R, there exists C > such that 

(4.7) \\kt',t)\{H'(M),H'-HM)) - ^' 

uniformly in t' and t. 

Proof. For u e 'io°°{M) we have u - Yjiei V^"- It thus suffices to take m/ e '^S'°°{AX), with supp(M,) c ft, for 
some / e J, and to prove that we have 

{d, + A(f'))f (,',o(M,) = it' - t)Lu,,){ui), Li,(,,,), e 'V\M), 
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and that Li^{t',t) satisfies (4.7) uniformly in t' and t. 

For concision we write q for q^u) here. Let us recall that qk is the Weyl symbol of Q{t) in the chart 
{®k,^k), k e K. We set p(j',){x,^) :- g^C "')9<''-*'f'. Making use of the assumption made on the chart 
(©<:(,), we consider the action of the change of variables k - o i/rj' on the operators ipj o 

p'J\,, Z).v) ° 'fj £ "^^{Sj) for i e J'l (see Figure 1). By Lemma B.2, we obtain 

with ^[°\) in *P''(A1) uniformly in t' and f. We then have 

(4.8) A(f') o = ° ( ^/ ° §" (f', X, D,) o o p;;: ,j(^, D,) o ) O (^,-('))*M; 

+ (r' - f)-R|,V,)«;, 

where ^['?,) is in *P^(A4) uniformly in f and f. 
From (4.2) we have 

d,'Pj.(i'.i)Ui = (A} ° (^> ° (qj{t, .)pjxi',i)y\x, £>x) ° 'it';) ° (tf^J^Tui, j G J",-, 
which we may write 

d,'Pjxt',t)Ui =(/'} q)(U X, Dx) o /?J(,,_,)(x, At) ° i^Pj) ° (^fYui 

where is in *P^(5)) uniformly in t' and f by Lemma 3.7. We choose € 'Tf^{6j) such that;^'^ is equal 

to one on supp(^y). We then have 

d,'Pj,(,',i)Ui ^ip*j ° {'fij ° q^j{t, X, D.,) oxj o /:'J(,,,,)(x, D.,) o^^. o ipj) o [ifrJ^Tui 
+ it' - 1)^0 j o r!^,,_^ o o ((Aj' )*«,-, e J-,-, 

recalling that q"'{t,x,Dx) is a differential operator, hence a local operator. Applying Lemma B.2 to Xj ° 

Pj(t'.tM' -^-v) o Afy we obtain 

where rJ^^, is in *I'^(A1) uniformly in t' and f and;fy = ^i/ry o ^l^'^^'j) Xi- Using again that q"(t,x,Djc) is a 
differential operator, we finally obtain 

(4.9) 5,P(,,„M,- = 2 ° ° ^' -^-^ ° ° ° ('J'mI))*"'' + (f' - tK'-»"'' 

where R[f,-^ is in *P^(A1) uniformly in t' and f. 

The operators in (4.8) and in (4.9) will contribute to the operator and we discard them 

(2) 

from the subsequent analysis. Observe that we may change the sums over j e J'i to sums over j e J'. 
(4.8) and in (4.9) since we only consider the action of the two operators on m,-. 
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Now that we have brought the analysis to the open set Qku), we shall consider and analyze the following 
symbol, o-(,'j), which corresponds to the operator Ci'^^^^Y ° (5, + A(f')) ° P(t',t) ° ^lio igiioring the operators 
and RfJ,^ as explained above: 



keeping in mind that we only consider the action of the associated operator on QP^^^.^yui whose support is 
compact and contained in ^k{i){0i)- We extend the symbol q{t, .) to R" x R" to obtain a symbol satisfying 
Assumption 2.1 like its counterpart in Section 2. We still denote by q{t, .) this extended symbol. We may 
then use global symbols in R". As in the proof of Proposition 3.6, we may replace q{t, .) by q{t' , .) by 
Assumption 4.2. For the symbol a-(f ,){x,^), this yields an error term of the form (r' - t)" A^^^ ^^{x,^), with 
^(f t)^^^^) in 5^ uniformly in t' and f, that will contribute to the operator Ln,ij). We thus discard this term 
in the subsequent analysis and we still denote by a-(f t){x,^) the modified symbol. 
We now set;^' = 1 - Zyej-'-' 'p] ™d we write 



Y (pjrq(t\)#''v]#'"Pi.',n#' 



which yields 

(4. 10) cr(,v)(x, ^) = -Y # Wc'.') ^'^i 

(4.11) + Y 

/, j e ■ ' 

From Weyl calculus [H6r85], since supp(;if) n suppiC^^^.^yuj) = 0, we find that the first term in the rh.s. of 
(4.10) can be written in the form (f' - t)A^^\x,^), with /l^'' in 5^ uniformly in t' and f, making use of the 
composition formula (1.6) and Lemma 2.3. 

Applying Proposition A. 1 (with k - 1), we find 

n 

^Y^ (2(5.>z)^^ - ^,5,„,(^j))%,p(,v, + (f - t)A^'-\x,^), 

m-1 
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with /l^^' in 5" uniformly in f and t arguing as in the proof of Proposition 2.4 in Appendix A (using 
Lemma 2.3 and Theorem 2.2.5 in [Kg81]). Therefore, we are now left with computing 

n 
n 

'"=1 l.jej^-' 

with /l*^* in 5^ uniformly in t' and f by the composition formula (1.6) and Lemma 2.3. Observing that 
the first term just obtained in fact vanishes, we finally have a-^,' ,){x,^) - {f - t)A{x,^) with A{x,^) in 5^ 
uniformly in t' and t. This concludes the proof. ■ 

4.4. Convergence and representation theorem. We observe that the energy estimate (3.5) also holds for 
the diff'erential operator A(f) on M (since the proof rehes on the Garding inequality which holds for positive 
elliptic operators on M). Combined with the (J?) stability result of Corollary 4.5 and the consistency 
estimate of Proposition 4.6, the energy estimate yields, as in Section 3, the following representation theorem 
through the convergence of 'Ws-^^t to U{t, 0), the solution operator of the parabolic Cauchy problem (4.4)- 
(4.5): U(t, 0) = limA,j,^o ^'^,t in the following sense. 

Theorem 4.7. Assume that A{t) satisfies Assumption 4.2. Then the approximation Ansatz "W^^j converges 
to the solution operator U{t, 0) of the Cauchy problem (4.4)-(4.5) in JXL^{M), i/^'^''(A1)) uniformly w.r.t. 
t as = maX()<j<Ar-i(f-'^' - ?*■'*) goes to with a convergence rate of order a{\ — r): 

- UitMkoM-'*') < CAli'-''\ f 6 [0, r], < r < L 

The operator ^i^t also converges to Uit, 0) in Lr{Q, T, X,{L^{Ai), L^{M))) with a convergence rate of order 
a: 

Furthermore 'TV-^.i,/ strongly converges to U{t, 0) in X.{L^{M), Lr{Ai)) uniformly w.r.t. t E [0, T]. 

Appendix A. Proofs of composition-like formulae 
We prove Proposition 1 . 1 and derive composition results for the symbol ph{x, f) - e 

A.l. Proof of Proposition 1.1. From Weyl Calculus we have 

(a rb){x, ^) = n-^" J ^'■^('•^•'•^•^>a(x + z, OKx + t, t) dz dt dr, 
® 



1 



22 HIROSHI ISOZAKI AND JEROME LE ROUSSEAU 

where S(z, ^, t, t, ^) = 2((r - ^, z> - <^ - ^, f>) (see [Hor85], p. 152). This yields 

(a #^'/7)(x, ^) - (ab)(x, ^) = n-^" J ei^(z,(AT,i) ^^(^ + + f, r) - a(x, OHx, r)) t/z t/^ c/f dr 

® 

(3) 

by a first-order Taylor formula. In the first (resp. second) term that we have obtained, we write 
Integration by parts w.rt. r and ( in the oscillatory integral yields 

n 2 ' 

— dtjOix + rz, ()d^.b{x + rf, r)) c/r c/z c/^ dt dr, 

2 J J e^^''^-'-^'^^o-(iD,,,D(),iDy,Dr))(a{x + rz,Ob(y+ dr dz d( dt dr _ , 

® 

which gives the result of Proposition 1.1 for k = O.To proceed further we integrate by parts w.rt. r and 
obtain 

(a r b)(x, f) - {ab){x, ^) = — — e'^<~'^'''^'^V((D^., D^), (D,„ D^)) OKy, t)) dz d( dt dr 

2 J ' \ I \y=jc 

" 2 ' 

+ Z ^ J(l - J e'^'^-^-'-'-MiD,, D^), (Dy, Dr)) {zjd,M(x + rz, Ob(y + rt, t) 

+tja(x + rz, Odx biy + rt, t)) dr dz d^ dt dr 

= !-cr{{D„D^),(Dy,D,)) [a{x,my,ri))\ 

+ ;r-2" ^1 j " / e'2(--.f.'.^.f) (cr((D„ D^), (Dy, Dr))f [a{x + rz, Ob(y + rt, t)) dr dz d(dtdT\_ , 

® 

which gives the result for k = I. Formula 1.6 then follows from induction by integration par parts w.r.t. r 
each time. ■ 

A. 2. From amplitudes to symbols. Here we give a formula of the form of (1.6) to compute the Weyl 
symbol of a tf/DO starting from an arbitrary amplitude. 

Proposition A.l. Let a{x,y,^) e S'"(R" x R" x R") be the amplitude of a iJ/DO A, i.e., 

Au(x) = (Ik)-" JJ e'<-^->''f>fl(jc, y, ^) u{y) dy d^. 
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The Weyl symbol b of A, i.e. A — (x, Dx), is then given by 



^ ;! \2 / \y=x 

1 ^ 

+ n-" J ^-^-jp- jj e^'^'-^-^^ i^^idx - dy, a{x + rz,y - rz, dr dz c/f . 



The proof is analogous to that of Proposition 1.1 given above. 

A. 3. Proof of Proposition 2.4. We prove the results for r/, #" /?/,. The results for p/, #">/, follow similarly. 

We first use Proposition 1 . 1 for k = 0: 

(A.l) 



- dx/h{x + rz, d^jPhix + rt, t) ) dr dz d( dt dr. 



@ 



By Lemma 2.3, we have 

with v/,(j) in 5 ^ and v*^* in 5 ' uniformly in h. We thus observe that the last term in (A. 1) can be written as a 

linear combination of terms of the form 

1 

® 

where vi_/, and vojt are respectively in 5'"'' and 5'"- uniformly in h with m\ +m2 — I + \ . Setting 

1 

Vh{x, 



we see that it is a multiple symbol in 5'"i>'"2(r2« x R" x R^" x R") and the term in (A.2) can be written as 



1 

\,x,y,y,^,ri) = J" vi,/,(rx + (1 - r)x,^)v2,h(,ry + {I - r)y,ri) dr. 



h J e'^<''^'''^'^VA(jc + z,x,y + t, y, t) dz d^ dt dr | 



® 

Applying Theorem 2.2.5 in [Kg81] twice (once for the integrations w.rt. z and t, a second time for the 
integrations w.r.t. t and (, recalling that S(z, t, t, ^) - 2((t - ^, z) - - ^, ?))) we obtain that the last term 
in (A.l) is of the form hA^i^\x,^), where /l*^''(x, ^) is in 5'^' uniformly in h. 
Similarly, by Lemma 2.3, we write 

dxjPh = h J Vh(j) , d^.ph = h 5 v^f , 

with Vh(j) in S ' and vj^"'* in S'^ uniformly in h. The same reasoning as above yields the last term in (A.l) is 
of the fovmh^Afix,^), where 

is in 5 uniformly in h. 
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To now treat the last equality in (2.2) we use Proposition 1.1 for^ - 1: 

(A.3) r/, rphix,^) = irhPh){x,^) + :^{r,„p,,}(x,^) 

2i 

^ 1 

- ld\,j:jh{.x + rz, 0d^.j^^ph(x + rt, t) + + rz, Od^.^^Phix + rf, t)) t/r c/z c/^ dt dr. 

Here, by Lemma 2.3, we write 

where v;,(j;i), vj^^^^^, and vj^'''*^' are respectively in 5 ^, 5 ' , and S ° uniformly in h. We can then conclude as above 
with Theorem 2.2.5 in [Kg81] and find the last term in (A.3) of the form hAf^ with if^ in S' uniformly in 
h. ■ 



Appendix B. Effect of a change of variables 

B. 1 . PseudodifFerential calculus results. We shall be interested in transformation formulae for Weyl sym- 
bols under a change of variables and apply them to the particular symbols we consider in Section 4. We let 
X and X be two open subsets of R" and let k : X ^ Xhea diffeomorphism. We shall study the eff^ect of the 
change of variables x i-> k{x) on the symbol x #" Ph #'''x the Weyl quantization, where ph - e^'^'', with q 
satisfying the assumptions made in Section 4 above and;^' e '^^^{X). 

We first consider general amplitudes before specializing to Weyl symbols. Let a{x, y, ^) be the amplitude 
in S"'{X X X X R") of A e *P'"(X) whose kernel is compactly supported. In particular, below, we shall 
consider a(x,3;,^) to be of the foTmx{x)x(y)a{x,y,^), with 5 e S'"{X xX xR"). With ( e ■^^^'"(R") equal 
to 1 in a neighborhood of we set 

aQix,y,^)^({x-y)a{x,y,^), and aUx,y,^) ^ {1 - ^{x - y)) a{x,y,^). 

If we set Ak - (k'^Y oAo k", then A,^ e ^'"'(X). In fact, for supp(^) sufliciently small, A^ - Ao,^ + Aoo,,^, with 
Ao^^it e ^''°°(X), and an amplitude of Ao,k is given by [GS94] 

(B.l) aoAx,y,0 = ao{K-\x),K-\y)J(iP'Hx,y)r'0\dsti'<-'yiy)\\detiF'ix,y)\-\ 

where K^^{x,y) - {K^i(x,y))i<i. i<„ is defined through 

n 

Af^'(x) - K^Uy) = ^ Ki^^(x,y)(x, - yi). 
1=1 

Note that k^^(x,x) - (K^^Yix) which implies that K^^{x,y) is indeed invertible in the support of oq when 
supp(^) is sufficiently small. Note also that 

(B.2) d,.F\x,y)\y^, = 5,. P'Hx,y)|,=, = ^d,,(K-'y{x), 
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Note that, for the operator Ac», we can regularize its kernel by integration by parts and use the amplitude 

n 

(B.3) at\x,y,^)^L''a^{x,y,^), with L = — 5— V - y,)5f„ A: 6 N, 

in place of aoo(_x, y, 

By Proposition A. 1 (with ^ = 1), the Weyl symbol of Ao,,f is given by 

(B.4) a,{x, ^) = e5«Ov,Of>-<o„Of»^j,_^(^^ ^)|^^^ ^ „-n jj e^'<-i-^> ao^x + z,x-z,0 dz 
= ao^„{x, x,0 + ^ (dx - dy, d^) a^x^ix, y, ^)|^_^^. 

=:a'».o(-t,f) =:a„,i(xS 
^ 1 

^ " (5)'/ U ^''^''^"^^(l - '^^^^^ - ^>'' ^f)'«o.^)(^ + '■^'3' - '•Z' * d^\^_^ . 


We now speciahze to an amphtude a(x, y, ^) given by the Weyl quantization, i.e., 

a(x,y,^)^X{x)x(y)b((x + y)l2S- 
To simplify some notation we set L = k K The symbol ff^.oCx, ^) is then given by 
(B.5) ff.,o(x, ^) = /7(L(x), f^'(L(x))^). 

Lemma B.l. The symbol ^) g/ven Z7j 

where 

n 

fk(x)^Yj^xMik(L(x))^ Yj {dl,x,'^,n){L{x)){d,^L,){x). 

l=\ l<m,/<n 

Proof. From the definition of a^ i in (B.4), and (B.l) we have 

aK,x(.xS = ^<5x-5,,5f) (Ar(L(^))A'(^3'))M(iW + L(3;))/2,^P'i(x,3;))-'a 

X |det(L)'(3;)||detP'i(x,y)ri)| , 

' \y=x 

where we have used that ( is equal to one in a neighborhood of the origin. From (B.2), we see that we 
need not take into account the spatial differentiations acting on the terms ^{K^^{x,y)Y^ . Similarly the spatial 
differentiations acting on the cut-off functions x{L{x)) and x{L{y)) cancel each other, and so do the spatial 
differentiations acting on the first variable of the symbol b. Note also that the absolute values for the last 
two terms can be removed before differentiation since their product yields 1 in the case y - x. To simplify 
the notation we set M - k^^{x, x). We thus obtain 

aK,i{x,^) = 2 (%^)(-L(x), ^K'{L{xm {M-% (5,. det(L'(x))) (detM)-'. 

l<j,k<n 
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From the multi-linearity of the determinant we find that 

(5,,det(L'(x)))(detM)-' = ^ ^;^/L(x)), 

l<p,l<n 

which yields 

fkix)^- 2 4(L(x)) (5,x;,(x)) ^;^(L(x)) = - 2 K',p(L(x)) {d„L'pj(x)) k'j,{L(x)) 

l<j,pj<ti l<j,p,l<n 

n 

^Yj^ML{x))), 

i=\ 

since k'{L{x)) L'{x) - Id^- ■ 

B.2. Application to tlie operator o p^{x,Djc) °X' We use the notation introduced above. In the case 
b = e"'"? = p,, then a(x, y, ^) - xix)x(y) e^'"'^^^^^^^^'^ is an amplitude for the operator A = x ° Phi^' ^.t) ° X 
with Weyl symbol a - x't^^Ph#^X- Making use of the form of the amplitude in (B.3), we see that 
Aco.K - hAoo.K with Aoo.K in ^'^{X) uniformly in h, using Lemma 2.3. 

We now focus on the operators Aq and Ao,k. From (B.5) and Lemma B. 1, the expression of the remainder 
term in (B.4) and using Lemma 2.3 we obtain 

n 

(B.6) a,{x, = xiLMfphiUx), ^k'{L{xM) (l - h'-^^ fk{x) V(L(x))^)] + ha,, 

with ff/f in 5° uniformly in h. Similarly, if we denote by the Weyl symbol of {iC^Y o q'\x, D J o a:*, we 
have 

n 

(B.7) q,{x, = qim, 'k'{L{xW + ^ ^ Mx) {d(,q){L{x), 'k'{L{x))0 + q., 

k=i 

where qi^ 6 5". We now prove that after the change of variables x i-» k(x), for the operator;^'opJJ (x, D^) ox = 
a^''(x, D v), we may use the symbol ;if(L(x)) #^^'x{L{x)) in place of Oi^ix, the pullback of a in the 

Weyl quantization, yet remaining within a first-order precision w.r.t. to the small parameter h. 

Lemma B.2. We set ph{x,^) = e-'"?'*-'-?). We have 

{((k-'Tx) r pH #" ((K-'rx)) (X, - = hAhix, 0, 
where Af, is in uniformly in h. 

Proof. We set v(x, ^) = i ^^^i fk{x){d^,q){L{x), 'k'{L{x))^. Making use of (B.7), we write 

Pkix,^) = /7„(L(x)//(L(x))^)e-'-<-^'«e-"^''-^'« 

1 

- Ph(L(x), 'k'(L(xW (l - hv(x, f) + (hvix, J e-''-(--f)(l - r) Jr) (1 + hiui{x, f)), 
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by two Taylor formulae, where fii is in 5^ uniformly in h. From Lemmata 4.1 and 2.3 we obtain that 



PhiLix), 'K'HixWUMx,^)- j e-'-'^'-^'^Hl -r)dr^ %2(^,^), 



with in 5'^ uniformly in h. From (B.6) we hence obtain 

with /i3 in 5° uniformly in li. We conclude the proof with the following lemma since pi, and pi, are of the 
same nature. ■ 

Lemma B.3. Let 4> e '(g';^{X). We then have 

(f>#'"ph#"'(p - 4>-ph = hAh, 

where Af, is in S uniformly in h. 

Proof. Since <p{x)(p{y)pi,{{x + y)/2,^) is an amplitude for the operator with Weyl symbol <p#'^'pi,#^'(p, by 
(B.4) we obtain 

(cprpi, r'4>){x, ^) = 7T-" jj e2'<-^.f-f>0(x - z) 4>{x + z) Phix, dz d^ 

1 

= (l>\x)ph{x, - \n-"Yj r rr^l - (-^C.^ + rz) 5?. - rz) 

- 2dxi4>{x + rz) d_^i^(f>(x - rz) + d\. .^^4>{x + rz) (f>{x - rz)) d^.^^j?h{x, () dr dz d(. 

We then conclude as in the proof of Proposition 2.4 in Appendix A by using Lemma 2.3 and Theorem 2.2.5 
in[Kg81]. ■ 
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